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ABSTRACT 



By means of direct N-body simulations and simplified numerical models, we study the formation and characteristics of the tidal tails 
around Palomar 5, along its orbit in the Milky Way potential. Unlike previous findings, we are able to reproduce the substructures 
observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities 
similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumu- 
lation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar 
clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, 
as previously suggested. Our models can reproduce the density contrast between the clumps and the surrounding tails found in the 
observed streams, without including any lumpiness in the dark halo, suggesting new upper limits on its granularity. 
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1. Introduction 



(Galaxy:) globular clusters: individual: Palomar 5 - Galaxy: evolution - Galaxy: kinematics and dynamics 



At a distance of 23.2 kpc from the Sun and of 16.6 kpc 
above the Galactic plane, Palomar 5 (Pal5) is a halo globular 
cluster (GC) characterized by very low total mass, Mtot = 5 
xlO^M©, a large core ra dius, = 24.3 pc, and low central 
concentration, c=0.66 (see lOdenkirchen et al.ll2002|). In the past 
decade, a number of observational studies (Odenkirchen et al." 
[200 1. 2003; Grillmair & Dionatos 2006; Jordi & Grebel 2010) 
have shown there are two massive tidal tails emanating from 
opposite side of the cluster, made of stars escaped from the 
system owing to the action of the Galactic tidal field. With 
an overall detected extension of 22° on the sky, corresponding 
to a projected spatial length of more than 10 kpc, Pal 5's tails 
are so elongated and massive that they contain more stellar 
mass than the one currently e stimated to be in the system 
itself (lOdenkirchen et al.l l2003l) . The resulting two tails are 
very narrow, with a proj ected, almost constant width of 120 pc 
(lOdenkirchen et al.l l2003i) . compa rable to the tidal radius of 
the cluster and kinematically cold (lOdenkirchen et al.ll2009l) . in 
agree ment with what can be expecte d from the disruption of 
a GC (ICapuzzo-Dolcetta et al. l l2005 |). Together with the tidal 
associated to NGC 5466 (iGrillmair & Johnson I l2006l) . Pal 5's 
tidal tails represent, so far, the most impressive example of the 
action of the Galactic tidal field on a GC. The large spatial 
extension of Pal 5's tails naturally permitted study of their fine 
structure. 

lOdenkirchen et al.l (l2003l) were the first to show evidence of 
inhomogeneities in Pal 5's streams, characterized by stellar 
density gaps (underdense regions) and clumps (overdense 
regions), particularly visible in the trailing stream, with the 
most evident and massive overdensities found between 100 and 
120 arcmin from the cluster center. The presence of density 
fluctuations along the stream has been successively confirmed 



by IGrillmair & Dionat"os ] (l2Q06h and iJordi & Grebel (l2QTQh . 
and the nature of these clumps has been long de bated. Already 
before the observations of Pal 5's tidal stream, ICombes et al.l 
(1999) showed the presence of small clumps in the tidal tails of 
simulated GCs, and attributed their formation to strong grav- 
itational shocks suflTered by the system at each disk crossing. 
When anal yzing a set of simulation s of GC s orbiti ng in a triaxial 
potential, [Capuzzo-Dolcetta et al.l ('2005) and Di Matte o~an 
(2005|), proposed an alternative explanation, suggesting for 
the first time that the clumps observed in their simulations 
were instead kinematic eflfects, related t o a local dece l eratio n 
in the motion o f stars along the tails. iKiipper et al.l (l2008l) . 
I Just e t al. (2009), and more recently Kii pper et al. r iisn2oT2h 
and ILane et al. (201Z) show that these kinematic eflTects and 
local decelerations correspond to "epicyclic cusps", where stars 
escaping from the cluster slow down in their epicyclic motion, 
as seen in the cluster reference frame. Because stars cross the 
tidal boundary with similar position and velocities, the location 
of such epicyclic cusps is similar for most escapers, generating 
visible overdensities. 

While all these works have successfully reproduced the presence 
of stellar substructures in tidal streams of GCs orbiting in an 
external potential, until now no simulation or model has been 
able to reproduce the only stellar clumps robustly observed 
in the tails of a galactic GC: those of Pal f]. In the most 
complete numerica l stud y of Pal 5's stellar streams realized 
so far, Dehnen al. (12001 have not succeeded in reproducing 
the clumpy nature of the stream, even after adopting a realistic 
Galactic model, and they propose that the observed clumps may 
be the eflTect of Galactic substructures not ac counted to i n their 
simulations, such as, dark matter subhalos. lYoon etal.1 (l201lh 



^ Clumps have been observed also alon g the GD-1 stellar str eam, 
whose progenitor was presumably a GC, see lKoposov et al.l (120101) . 
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have recently discussed, by means of simplified numerical mod- 
els, the absence of clumpy substructures in Pal 5's tails when 
the cluster is evolved in a smooth dark matter halo, instead 
showing that stellar overdensities and gaps between them can 
be generated if the gravitationa l eff^ect of dark s ubhalos on the 
streams is taken into account. ICarlberd (l20T2h accounts for 
the possibility that epicyclic motions could explain the density 
variations observed in the inner part of the tails, but excluded 
the possibility they could reproduce those beyond about 1 kpc 
(2.5°) from the cluster center. From the density of underdense 
regions in observed narrow streams (Pal 5, among others), 
he estimates a total halo population of about 10^ dark matter 
subhalos with masses greater than lO^M©. 
The results we present in this Letter reproduce, for the first 
time, the inhomogeneous properties of Pal 5's streams without 
invoking any lumpy dark halo. As described in the following 
sections, epicyclic pileup of stars can be responsible for gen- 
erating these substructures, even at significant distances from 
the cluster. We also exclude either gravitational instabilities 
or episodic mass loss from the cluster as processes that may 
explain their formation. 



2. Models and initial conditions 

2.1. N-body simulations and restricted three-body methods 

We have studied the formation and characteristics of Pal 5's tidal 
tails by means of N-body, direct summation simulations, and 
simplified numerical models. To run the simulations we used a 
modified version of NBSymple (ICapuzzo-Dolcetta et al 
a high performance direct A/^-body code implemented on a hy- 
brid CPU-hGPU platform by means of a double-parallelization 
on CPUs and on the hosted graphic processing units (GPUs). 
The eff'ect of the external galactic field is taken into account by 
using an analytical representation of its gravitational potential 
(see ^2.21) . We ran several simulations, varying, among other 
parameters, the value of the softening length adopted to smooth 
the mutual interactions. In this Letter we present the results of a 
simulation with a softening length s = 0.03 pc, but we want to 
emphasize that in all these simulations, the presence of substruc- 
tures was detected. The time integration of the particles trajec- 
tories was done using the second-order leapfrog method, with a 
time step = ^s^/(GM), where M is the total mass of the clus- 
ter. To fully exploit the GPUs power the force was computed in 
emulated double precisionWith these choices, the average rela- 
tive error per time step was less than 10"^^ for the cluster both 
isolated or in the Galactic potential. 

The result of N-body simulations were compared to simpli- 
fied, restricted three-body models, where we used the same ini- 
tial conditions adopted to run the simulation. However, at each 
time step, instead of evaluating the mutual interactions among 
all the N particles in the system, we evaluated the interaction 
with the global gravitational potential of the cluster and with the 
Galactic potential for each particle. This approach allows the 
main characteristics of stellar streams to be reproduced (orien- 
tation, spatial distribution and extension, presence of overdense 
and underdense regions) without having to make CPU-time con- 
suming N-body simulations. It represents, of course, a simplified 
model that does not take the evolution of the internal parame- 
ters of the system (its mass and concentration) into account and 
thus must be taken as a first-order approach in the study of tidal 
streams. 
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Fig. 1. Map of the surface density of stars of the simulated cluster 
in the plane of the sky. The solid red line represents the local 
orbit. 

Finally, similar to lLane et al.l (l2012h , we produced streak- 
lines that allow for a simple visualization of the spatial distri- 
bution, at the present time, of stars that have escaped from Pal 
5 in the last Gyrs. To do so, we integrated the orbit of a point 
mass representing Pal 5 bary center, and, at each time step, re- 
leased two particles at distance r^ax = 120 pc from it, along the 
instantaneous Galactic center and anticenter directions, and with 
a velocity equal to that of the cluster. The motion of these esca- 
pers is integrated in time, taking both the eff'ects of the GC and 
Galactic potential into account. 

2.2. Globular cluster and Galaxy models 

The Galactic potential is from I Allen & SantillanI (Il99ll) . whose 
model consists of a three-component syst em: a spherical cen- 
tral bulge and a flattened disk, both in the iMiyamoto & Nagail 
(1975) form, plus a massive smooth spherical halo. We did not 
include a barred potential, since its eff'ect on the orbit of Pal 5 
has been shown to be negligible (I Allen et al.ll2006h . 

For th e GC we adopted t he parameters of the "model A" de- 
scribed bv lDehnen all (120041) . because, according to the authors, 
this model is the one that best fits the stellar num ber density 
of Pa l 5 at the end of their simulations (see Fi g. 12.jDehnen al.l 
|2004|) . It consists of a single-mass King model (lKingll966h with 
tidal radius = 56 pc, Wq = 2.75, and Mq = 2 x 10"^ M©. To 
generate positions and velocities of th e N = 15360 p articles that 
represent the cluster we used NEMO (iTeubenll 19951) . 

For the orbital initial conditions, we first integrated a test 
particle with Pal 5's present position and v elocity back \ yard i n 
time for the same time interval as chosen by lDehnen al.l (|2004|) . 
i.e. 2.95 Gyr, in the Galactic potential, and obtained the ini- 
tial conditions of the cluster barycenter. Then we translated the 
positions and velocities of the single-mass King model previ- 
ously described to coincide with those of the cluster barycenter 
2.95 Gyr ago, and we integrated the whole cluster 2.95 Gyr for- 
wa rd in time. The current pos ition and velocity of Pal 5 are given 
by lOdenkirchen et aP (l2003l) . In particular the position of Pal 5 
in the Galaxy is ix,y,z) = (8.2,0.2, 16.6) kpc; here x, y, and z 
denote righthanded Galactocentric Cartesian coordinates, where 
the Sun has coordinates (-8.0, 0.0, 0.0). We assumed a radial ve- 
locity Vr = -44.3 km/s (observer at rest at the present location 
of the Sun), a tangential velocity (as seen by observer at rest at 
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Fig. 2. Upper panel: The trailing tail at the end of the N-body 
simulation. The contour lines in the tails refer to 740, 2200, 
4800, 6300 and 7400 Mo/deg^. Bottom panel: Linear density 
of the trailing tail shown in the upper panel, as a function of the 
right ascension a. The cluster is on the right of the panel. Error 
bars are Poisson errors. 



the present location of the Sun) Vt = 90 km/s with a P.A.= 231° 
with respect to the direction pointing to the northern equatorial 
pole, and 280° with respect to Galactic north as implied by the 
fit of the local orbit (lOdenkirchen et al.ll2003l) . These data yield 
a spatial velocity with components (-40.7, -89.3, -21.0) km/s 
in the Galactic reference frame mentioned above. 



3. Results and discussion 

After 2.95 Gyr of evolution, 15 disk crossings and 11 pericenter 
passages, which every ~ 270 Myr bring the cluster to distances 
of ~ 5.9 kpc from the Galaxy center, the N-body simulation 
shows that the cluster has lost most of its initial mass (88%), 
now redistributed into two long and narrow tidal tails. The ap- 
parently smoothness of the tails hides a complex and inhomoge- 
neous stellar distribution that becomes evident only when plot- 
ting the isodensity contours, as in Fig. [T] for the whole stellar 
stream. This plot reveals that indeed the tails are characterized 
by several substructures, and it shows that currently the dens- 
est are all localized in the trailing tail at distances between 0.5° 
and 3.5° from the cluster center. In this portion of the tail, our 
N-body simulation predicts there are two prominent density en- 
hancements (see top panel of Fig. O: the first located very close 
to the cluster (a < 230°) and the second starting at a=230.5° 
and ending at a=232°. The second region shows an underdense 
region in it, whose extension is ~ 0.5°. All these characteris- 
tics, at the same location, are also found in the observed streams 
(cf, for example. Fig. 3, in ^denkirchen et al. 2003). Not only 
does the position of these clumps closely resemble the loca- 
tion of those observed in Pal 5's tails, but their linear densities 
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Fig. 3. The trailing tail of the cluster in the restricted three-body 
model. The contour lines in the tails refer to 470, 940, 1410, 
1970, 2500 and 3300 Mo/degl 



(bottom panel of Fig. [2]) - 2-3 times above the density of the 
surrounding streams - and the presence of underdense regions, 
with sizes of 0.5°- 1°, also agree with those measured for Pal 
5 (see Fig. 4. lOdenkirchen et aLll2003l) . These substructures are 
also visible in the isodensity contours of the simplified numer- 
ical model (in Fig. [3]), at few degrees from the cluster center. 
When compared to the full N-body simulation, the diff'erences in 
the exact location and amplitude of the clumps in the restricted 
three-body problem occur because the former does not take the 
internal evolution of the cluster into account. The gravitational 
potential of the cluster is kept fixed with time, and thus the in- 
stantaneous mass loss that aff'ects the density of the innermost 
regions of the tails can only be captured approximatively un- 
like in the N-body simulation, where these phenomena are mod- 
eled self-consistently. This ultimately leads to some diff'erences 
in the density distribution along the streams, with clumps arising 
in the simplified numerical model that are less massive and less 
dense. Nevertheless, both approaches unequivocally show sub- 
structures in models of Pal 5's tidal tails, a few degrees from the 
cluster center (see Appendix [A] for the dependency of the sub- 
structures on the orbital parameters). 

These regions of over- and underdensities are due to the 
epicyclic motion of star escaping from the cluster: stars lost 
at diff'erent times redistribute along the tail following a com- 
plex path, as shown by the streaklines in Fig. |4l This kine- 
matic process, described and studied in a number of papers 
(KiiDper et al.""2008i lJust et al.ll2009l: iKiipper et al.ll20 Id 120121: 
Lane et al. 2012), is applied here for the first time to reproduce 
observed stellar streams, including both the position and the in- 
tensity of the observed stellar inhomogeneities. Also the absence 
of overdensities in the trailing tail at distances greater than 4°-5° 
is simply because we have integrated Pal 5's orbit for the last ~3 
Gyr. With this choice, at the current cluster position, most of the 
escaped stars do not reach projected distances over 4°-5° from 
the center of the cluster 0. If the orbit is integrated over much 
longer times (~8 Gyr), streaklines show that epicycle loops also 
form at distances of several kpc from the center. This process is 
thus able to form over- and underdensities in several regions of 
the tails, at a variety of distances from Pal 5. 



^ The length of tidal tails depends on the location of the GC along 
its orbit: for example, at the last pericenter passage, 160Myr ago, our 
models show that stars were spread along a 12 kpc stream, 2.2 times 
longer than the extension of the simulated Pal 5's stream at the current 
position. 
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Fig. 4. Streaklines which show the positions, at the present time, 
of stars that have escaped from Pal 5 in the last 2.95 Gyrs. 
Different colors correspond to stars lost between 0. and 0.5 Gyr 
ago (red dots), between 0.5 and 1.0 Gyr ago (green dots), be- 
tween 1.0 and 1.5 Gyr ago (blue dots), between 1.5 and 2.0 Gyr 
ago (magenta dots), between 2.0 and 2.5 Gyr ago (cyan dots), 
between 2.5 Gyr and 2.95 Gyr ago (black dots). 



Moreover, that these substructures are also reproduced in 
models where the mutual gravity is neglected rules out Jeans 
instability as a possibl e mechanism responsib l e for th eir forma- 
tion, as suggested by lOuillen & Comparettal (|2010|). This hv- 
pothes is has been already questioned by 'Schneider & Moore 
(1201 Ih . and here we confirm that it can be excluded as the mech- 
anism at the origin of the substructures observed in Pal 5. 

The association of overdense regions along the tails to the 
location of epicycles (upper panel of Fig.[2l Figs[3l and |4]) also 
rules out the possibility that they are due to an episodic mass loss 
from the cluster (see Appendix |B] for further details). 



4. Conclusions 

We have studied the formation and characteristics of the tidal 
tails of the GC Pal 5 along its orbit in a smooth Milky Way po- 
tential by means of N-body simulations and simplified numerical 
models. Our main results can be summarized as follows: 

- The trailing tail of the cluster is more inhomogeneous than 
the leading tail. 

- The densest overdensities are found in the trailing tail at dis- 
tances between 0.5° and 3.5° from the cluster center. 

- In this portion of the tail, our N-body simulation predicts 
there are two prominent density enhancements, the first lo- 
cated very close to the cluster (a < 230°) and the second 
starting at a=230.5° and ending at a=232°. The second en- 
hancement shows an underdense region, whose extension is 
~ 0.5°. These characteristics, at the same location, are found 
in the observed streams. 

- The linear density in the clumps is two to three times above 
the density of the surrounding streams and the underdense 
regions in between clumps have sizes of 0.5°- 1°, comparable 
to those measured for Pal 5. 

- Pal 5's clumpy streams can also be fo rmed in a smooth 
dark halo, unlik e in previous findings (lOuinn et al. I l2008t 
lYoon etaDl201lh . as a simple consequence of the epicyclic 
motion of stars in the tails. 

- Jeans instability and episodic mass loss can be ruled out as 
mechanisms at the origin of the observed substructures. 



To conclude, our models can reproduce the density contrast 
between the clumps and the surrounding tails found in the ob- 
served streams, without including any lumpiness in the dark 
matter halo. At the same time, many works (among others, see 
llbata et al.ll2002H Johnston et al.l2()02h have shown that the struc - 
ture, as well as the kinematics dSiegal-Gaskins & Vallurill2008l) . 
of stellar streams can be very sensitive to heating by encounters 
with massive dark-matter subhalos. It would be thus interesting 
to reconsider the impact of these substructures on the character- 
istics of Pal 5's stellar streams in light of these new results and, 
in particular, to set new upper limits on the granularity of the 
Milky Way dark halo. Finally, future models may want to take 
also inhomogeneities in the Milky Way stellar disk into account 
(molecular clouds, spiral arms) to derive a complete picture of 
the morphological evolution of stellar streams. 
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Appendix A: Dependency on the orbital parameters 

In Section [3l we have seen that epicyclic motion of stars in the 
tails generate overdensities mostly in the trailing tail of Pal 5, at 
distances of few degrees from the cluster center. Here we show 
how sensitive this finding is to a change in the orbital parameters. 
To this aim, we have run several N-body simulations changing 
both the radial, v^, and tangential, Vt velocities with respect to 
those of the standard orbit (we recall the reader that for the stan- 
dard orbit, Vr = -44.3 km/s and Vt = 90 km/s). In particular, for 
the tangential velocity, we adopted the extreme values, 80 km/s 
and 110 km/s, beyond which the orbit fails to give a good rep- 
resent ation of the observations, as shown by lOdenkirchen et aH 
(l2003h . and also two intermediate values similar to those of the 
standard orbit, = 89 km/s and = 91 km/s. For the radial 
velocity we arbitrarily varied the best fit parameter by ±2 km/s, 
since the error on this value is smaller than th e one aff'ecting 
the tangential velocity (lOdenkirchen et al.ll2003l) . In all the mod- 
els presented in this Section, the position angle P.A. has been 
kept fixed to = 231° and the internal parameters of the clus- 
ter have not been changed. In Figs. lA.ll and IA.2l the results of 
this analysis are presented. The orbits with Vr = -44.3 km/s 
and Vt = 110 km/s or Vt = 80 km/s are significantly diff'erent 
from the standard orbit. On the first orbit (which has a pericenter 
-1.2 kpc greater than that of the standard orbit) the cluster loses 
a small percentage of its mass, and the tails are characterized by 
a low density and a nearly flat stellar distribution, with density 
variations within the error bars. On the other hand, although the 
cluster with Vr = -44.3 km/s and Vt = 80 km/s is destroyed by 
the tidal interaction with the Galactic field (the pericenter of the 
orbit is 4.7 kpc, a factor -1.3 smaller than that of the standard 
orbit), clumps are clearly visible between 229° and 231°. Thus, 
for variation of the tangential velocity of the order of 10 km/s the 
density of the tails changes significantly. For smaller changes of 
the tangential and radial velocities the mass loss rate is closer 
to that of the standard orbit and stellar overdensities are present 
in the trailing tail, with similar positions and amplitudes. This 
also confirms that the overdensities found for the standard orbit 
are not due to random fluctuations, but rather robust characteris- 
tics of this portion of the tail, since for small changes around the 
standard values they are still present, and generate a profile that 
closely resembles that of the standard orbit. 

Appendix B: Why are clumps not caused by an 
episodic mass loss? 

In their modeling of Pal 5's tidal tails, iDehnen aL I (l200l have 
already shown that at any given distance from the cluster cen- 
ter, stars lost from the cluster at diff'erent times can be found. 
We confirm this finding in Fig. lA.3l where we plot the time stars 
escape from the tidal boundary 0as a function of the right ascen- 
sion a, which describes the spatial extension of the streams. At 
any given value of a corresponds many diff'erent values of tioss' 
this is especially the case for particles in the inner part of the 
tails, at distances of 3° - 4° from the cluster center. Moreover, the 
distribution in the plane a - tioss resembles a tree, with branches 
and trunk corresponding, respectively, to stars lost at each peri- 
center passage (plotted in Fig. IA.3I with diff'erent colors), and 
in between two consecutive pericenter passages. Due to their 
higher velocity dispersion when they escape the cluster, stars lost 
at the pericenter passage are rapidly spread over a larger exten- 
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Fig.A.l. Upper Panel: linear density of the trailing tail, as a 
function of the distance from the cluster center, for diff'erent or- 
bits having the same radial velocity of the standard orbit (S.O., 
black line), i.e. -44.3 km/s, and Vt = 80 km/s (blue line), 
= 89 km/s (magenta line), Vt = 91 km/s (green line), or 
Vt = 110 km/s (red line). Lower Panel: The same but for or- 
bits having the same tangential velocity of the standard orbit 
(S.O., black line), i.e. 90 km/s, and Vr = -42 km/s (blue line) 
or Vr = -46 km/s (red line). The typical Poisson's uncertainty is 
shown. 

sion of the tidal tails than stars lost in between consecutive peri- 
center passages, which are systematically, for every tioss, closer 
to the cluster center than stars lost at pericenters. This behavior 
in the redistribution of stars lost in diff'erent phases of the clus- 
ter orbit clearly demonstrates that any temporary accumulation 
outside the tidal boundary of stars lost at the pericenter passage 
is rapidly cleaned out and thus cannot be a mechanism able to 
produce stellar overdensities in tidal streams at several degrees 
of distance from the cluster center. 



^ For simplicity we adopted the current tidal radius of the cluster as 
reference value. 
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Fig.A.2. Isodensity contour plots for the trailing tail of the cluster for different N-body simulations having orbits with Vr = 
-44.3 km/s and Vt = 80 km/s (upper left panel), Vr = -44.3 km/s and Vt = 89 km/s (upper right panel), Vr = -44.3 km/s and 
Vt = 91 km/s (middle left panel), Vr = -44.3 km/s and = 1 10 km/s (middle right panel), Vr = -42 km/s and Vt = 90 km/s (bottom 
left panel), Vr = -46 km/s and Vt = 90 km/s (bottom right panel). See upper panel of Fig. [2l for a comparison with the isodensity 
contours of the standard orbit 
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Fig. A.3. Time, tioss, at which escaped stars crossed the tidal boundary of the cluster as a function of the right ascension a of stars 
in Pal 5's tails. The right ascension has been evaluated at the current epoch (see Fig. [T] for comparison), tioss is zero for stars lost 
from the GC at the current epoch, and attends the lowest values for stars escaped from the GC at the beginning of the simulation. 
Each color in the plot corresponds to stars lost at different pericenter passages, while grey points correspond to stars lost between 
two consecutive pericenter passages. 



